STUDY OF THE VACUUM MATRIX ELEMENT 
OF PRODUCTS OF PARAFIELDS 



O.W. Greenber^ 

Center for Fundamental Physics 
Department of Physics 
University of Maryland 
College Park, MD 20742-4111 
and 
A.K. MishreE] 
Institute of Mathematical Sciences 
Chennai 600 113, India 
University of Maryland Preprint PP-09-049 

[^Xiv:n9n9.4n69l 

Abstract 

We study the vacuum matrix elements of products of parafields using 
graphical and combinatorial methods. 

1 Introduction 

We can calculate matrix elements of parafields using either the trilinear commutation 
relations of Green, or using the Green ansatz, also due to Green in the same paper [1] . 
We recently gave a path integral quantization of parastatistics [2], primarily to 
show that parastatistics can be quantized using path integrals. We found that, in 
both the canonical and path integral formalisms, the Green ansatz is convenient 
for calculating quantities in parafield theories. This is in contrast with the trilinear 
commutation relations that are cumbersome to use for calculations. 
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As far as we know, there is no general calculation of the vacuum matrix element 
of an arbitrary product of parafields in the literature. Such a calculation will be 
useful in applications of parastatistics, such as a calculation of the partition function 
and its connection with the counting of partitions in number theory [3] . In addition 
this calculation presents combinatorial problems that are interesting in their own 
right. In this paper we present the explicit calculation of the vacuum matrix element 
of products of < 6 scalar parabose fields (and spinor parafermi fields) of any order 
p. We also give examples of terms in the vacuum matrix elements of = 8 and 
A^ = 12 parafields. In addition we give methods that are helpful to calculate the 
vacuum matrix elements to arbitrary order. 

2 The vacuum matrix element of a product of N 
parabose fields of order p 

The matrix element of an arbitrary product of free bose fields is the sum over all 
products of contractions, each of which is a two-point function with coefficient one. 
This is the "factor pairing theorem." [H [S] 

For the corresponding case with parabose fields, not all coefficients are one; 
indeed the main difficulty in giving a general result for parabose fields is to calculate 
the coefficient of each term. These coefficients are positive or negative integers or 
zero. 

When we work in x-space we use the Pauli- Jordan commutator function which 
is odd in x, 




(1) 



and 




(2) 



where 





For completeness we also define the even function of x 




(4) 
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The analogous formulas for free spinor parafermi fields are 

iS{x) = {if^ + m)iA{x), (5) 

tS{x) = S'^+\x) - S^-\x) (6) 

where 

S(^\x) = (iy, + m)iA^^\x). (7) 

S^^\x) = S'^+^x) + S'^~\x). (8) 

The formula that we gave in our previous paper, 

0(")(xi)0(^)(a;2) = {25^13 - l)<P^^\x2)<P^''\x,) + 5^ptA{x, - x^) (9) 

is very useful for this calculation. In the present context, 

(25,,,^, - 1)2 = 45„,„^. - 45,,,^, + 1 = 1. (10) 

This is obvious without calculation, since 

26a,a, - 1 = 1, i = J, and - 1, 2 7^ j; (11) 

i.e., this expression has the value ±1 and thus its square is 1. This simple identity 
plays an important role. It guarantees that the result is independent of the details of 
the procedure used to evaluate the matrix element and also is useful in evaluating the 
sums over the a's. (Another way to write this expression is 25^^.^^^- — 1 = (—1)''"''+^.) 
Here the {(j)^"^} fields are the Green components associated with the parabose fields 
{$}. The analogous formula for parafermi fields is 

^'("Ha^O^^^Ha^a) = (1 - 25,^)^(^)(x2)V'^"Ha;i) + M5(xi - x^) (12) 

We choose the following systematic way to evaluate a vacuum matrix element, which 
we illustrate for the four-point function]^ 

(0|$(X1)<I>(X2)$(X3)$(X4)|0) = 

(0 (xi)0("^) (X2)0("^) (^3)0^"^^ (a;4) |0) (13) 



■^For completeness we record the result for the two-point functions, 

< 0|$(xi)$(a;2)|0 >= pAW(x2 - xi) 



< 0|^'(xi)i'(x2)|0 >= pS'^+\x2 - Xi) 
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Let ( = + act to the left on the vacuum, so that the positive 

frequency part, i.e. the creation part (f)^"^^~\ annihilates the vacuum on the left. 
Then in Eq. (1) the iA function will be replaced by the A*^"*"^ function. Now move 
i^i^Q right until it annihilates the vacuum on the right. This field will 
contract with each of the three fields on its right, giving a Saia, factor multiplying 
a matrix element with a product of two fields plus a term in which has moved 
one transposition to the right multiplied by a (26 — 1) factor. The remaining two- 
point functions can now be evaluated in a similar manner, except that the second 
type of term just described with the [26 — 1) factor will be absent. Where a factor 



C^SaiOj — 1) occurs, we replace it by 1 according to Eq.(lO). For our illustrative 
case of the four-point function, this procedure, which we call the first step in our 
general discussion below, yields 

(0|$(X1)$(X2)$(X3)$(X4)|0) = 

5aia2A'^^\xi - X2)5a3a4A^^\x3 - X4) + 

a's 

(2(5„ia2 - l)5ai„3A(+)(xi - X3)6a2aiA'-^\x2 - X4) + 

{2Saia2 - l)(25„i„3 - l)5aia4A^^\xi - X4)Sa2a3^^^\x2 - X3). (14) 

For the second step of the calculation of this four-point function, each term has two 
different factors 6a,aj that we use to do two of the sums over the a's. By renaming 
the dummy a's, if necessary, we reduce the result to a sum over the first two a's. 
The result is 

(0|$(xi)$(x2)$(a;3)$(x4)|0) = 

J2 [^^""^ i^l - 2:2) A W (X3 - X4) + 

(26^,^, - l)A(+)(xi - X3)A(+)(x2 - 0:4) + 

(25„,„, - 1)(25„,„, - 1)AW(xi - X4)A(+)(a;2 - 0:3)], (15) 



where we used the identity Eq.(lO) in the third term. In the third and final step. 
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we do the sums over the two remaining a's. This gives the result 



(0|$(xi)<l>(a;2)$(x3)<f(x4)|0) = 

p2^(+) (a;^ _ a;2) a(+) (xs - X4) + p{2 - p)A^^^ {xi - X3) A^^^ {x2 - X4) + 
p2AW(xi - X4)A^+\x2 - xs). (16) 

For charged spinor parafermi fields the calculation is similar, except we keep 
only contractions between and field ip and its conjugate ■0. 

In the general case the first step of this procedure reduces the matrix element 
with N — 2n fields to a sum oi N — 1 terms, each of which will have a product 
oi N — 2 fields. The entire procedure yields {N — 1)!! terms, each of which has 
the product of n two-point functions A*^+) multiplied by the product of n 5's and a 
product of (25 — 1) factors. These terms will be summed over the N a''s from 1 to 
p for each a. 

We represent each term by a graph containing N points and n directed lines 
or "links" below the points representing the contractions. We can arrive at these 
graphs directly by forming all possible {N — 1)!! distinct combinations of n-pairs 
involving N — 2n integers 1 to N. In the graphs, these integers are represented by 
N points, numbered in ascending order. Subsequently, we connect the pair of points 
representing the integers in a pair, by a line in the graph. As each combination has 
n pairs, n such lines are required. We choose the line (or link) to go from the lower 
numbered point to a higher numbered one. 

Next we provide an algebraic correspondence for these graphs. We start with 
a given graph and associate a factor 5aiajA^^\xai — Xa^) with the line linking the 
points i and j . We repeat this process for each of the n links, and form the product 
of all factors Saf,aiA^^\xa^ — Xai). This resulting term has to be multiplied by a 
product of {(2(5 — 1)} factors. The number of crossings of the links corresponds to the 
number of {25 — 1) factors in the product. (All of the links must lie below (or above) 
the line containing the N points.) To arrive at their specific values, start with the 
line originating from the point 1, and terminating, say, at point t. If this line crosses 
other lines originating from points ^1,^2, •••^g, all lying to the left of the point t, we 
get a product of {{25 — 1)} factors, namely, {{25aiaj — l),^i < j < tq}- If the line 
originating from point 1 does not cross any other line, then no {25a■^aj — 1) factor 
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is introduced. Repeat this process for the hnes originating from points 2,3, ...,n, 
and for each new crossing encountered in this process, introduce a {25 — 1) factor. 
Finally, take the product of all these {{26 — 1)} factors, multiply this quantity by 
the product of n-pairs of ^^.a^ — Xa^) factors obtained earlier, and sum the 

resulting expression from 1 to p over ai,a2, ■■■,ckn- 

As an illustration, we consider here the N = 6 case. Generating all the five 
0-crossing, six 1-crossing, three 2-crossings and the single 3-crossings graphs for N 
=6, and substituting the associated coefficients derived in the earlier paragraph, the 
six-point function is 

(0|$(X1)$(X2)$(X3)$(X4)$(X5)$(X6)|0) 

= p^{A(+)(a;i - X2)A^^\x3 - X4)A^^\x5 - xq) + A(+)(a;i - X2)A^^\x3 - xq)A^+\x4 - a 
+A(+)(xi - X4)A(+)(a;2 - a;3)A(+)(a;5 - xq) + A(+)(xi - Xe)A^^\x2 - X3)A^^\x4 - x^ 

+AW(xi - X6)AW(x2 - X5)A^+\X3 - X4)} + 

p\2-p) 

X {A(+)(xi - xe)A^+\x2 - X4)AW(x3 - x,) + A^+\x^ - X5)A^+\x2 - xs)A^+\x4 - x^) 
+A(+)(xi - a;3)A(+)(a;2 - a;6)A(+)(a;4 - x^) + A(+)(a;i - a;3)A(+^(x2 - a;4)A(+)(x5 - Xe 
+A(+^(xi - a;2)A(+)(a;3 - X5)A(+)(x4 - xq) + A'^+\xi - X5)A^+\x2 - X6)A^+\x3 - X4 
+p{2-pf 

X {A(+)(xi - X5)AW(x2 - X4)AW(x3 - xe) + A(+)(xi - X4)AW(x2 - X^)A^+\X3 - X5) 
+A(+)(xi - X3)A^^\X2 - X5)AW(x4 - Xe)} 
-p{p^ - 6p + 4){A(+)(a;i - 0:4) A^ (2:2 - a;5)A(+)(a;3 - xq)] 

(1 

Next, we give the following 4 cases corresponding to all possible 5-crossings 
graphs for an 8-point function as illustrations. The algebraic expressions for graph 
(1), graph (2), graph (3) and graph (4), drawn in figures (1-4), are 

{25a-,,a2 - l)(25aia3 - 1) (2(5^104 - 1) (25^203 " I)(25a2a4 " 1) X 

a\,a2,...,as 

(5c,i„5A(+)(a;i - x^)5a,,a,^A''''\x2 - xq) x 

5o,^,aA^+\x3 - X8)5a4a,A(+)(x4 - X7), (18) 
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^ ] [(2(5aj^a2 -'-)(25aia3 ^){'^^a2,a3 
A(+)(x3-X7)5a,a3A(+)(x5-X8)], 



1)(25, 



1)(25, 



1)X 



(19) 



[(25aia2 - l)(25a^a3 " l)(25aia4 " I)(25a2a4 " I)(25a3a4 " 1) X 



ai,Q;2r--)0:8 



A^+)(a;i - X5)5a2aT^^^\x2 - xr) x 

A(+)(x3-X6)5«4asA(+)(x4-X8)], 



(20) 



and 



[(25aia3 - l)(25aia4 " I)(25a2a3 " I)(25a2a4 - I)(25a3a4 " 1) X 

ai,a2,...,as 



Saiae^^'^Hxi - Xe)Sa^a^A''+\x2 - X5) X 
A(+)(X3 - X7)5^,as^^+\X4 - Xs)], 



(21) 



respectively. 

Since there are n (5-f unctions in each one of the {N — 1)11 terms, in the second 
step we reduce the expressions to sums over n a's. We now define a graph with n 
vertices and as many hnks between vertices as there are crossings in the diagrams 
with N points discussed above. A hnk connects a pair of vertices if fines starting on 
eacfi of tfie pair intersect. Tfie most "saturated" grapfi witfi tfie maximum number, 
n{n — l)/2 finks corresponds to a simplex in n-dimensional space. This also implies 
that the coefficients contain at most a product of [n{n — l)/2] {26 — 1) factors. 
In the final step wc do the sums over the remaining n a's. The total number of 
summations to be carried out is (A^ — 1)!!. Obviously for a graph with m crossings 
(1 < m < n(n — l)/2), the summand is a product of m {26 — 1) terms. When m 
= 0, summand is unity, and the cocflicicnt equals p". For m = 1, the summand is 
(25aj,aj ~ 1)) and the associated coefficient takes the value p"'^'^{2 — p). When m = 
2, all possible values of the summand, through relabeling the dummy indices {ai}, 
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get reduced to two forms, viz, {2Sa-,,a2 - l)(25ai,a3 - 1), and {2Sa-,,a2 - I)(2(^a3,a4 - 1)- 
The coefficients corresponding to both the summands are p"~^(2 — p)^. For m 
— 3, N — 6, only one graph with three crossings exists; having the summand as 
(25ai,a2 ~l)(2(^ai,a3 — I)(25a2,a3 ^nd the Coefficient is — — 6p+4). In general, 
when m takes the maximum value n{n - l)/2, the is ni<j;i<i<n-i;2<j<n(2<^ai,a^ - 1), 
and there is only one graph corresponding to this most saturated case. 

As a step towards finding a general method to evaluate the coefficients, we can 
simplify the problem of (A^ — 1)!! summations for N number of fields by enumerating 
the number of graphs Tjvm having m numbers of crossings. As noted earlier, < 
m < n{n — l)/2. In the following table, these values are listed for 2 < A^ < 10. The 
* in a box implies that the m has reached its maximum permissible value. 
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TABLE I 



m 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


/AT" 1 N 1 1 

= (A* — 1)!! 


■J- 2m 


1 * 
i 






















i 


Tim 


2 


1* 




















3 


T 
6m 


5 


6 


3 


1* 
















15 




14 


28 


28 


20 


10 


4 


1* 










105 


TlOm 


42 


120 


180 


195 


165 


117 


70 


35 


15 


5 


1* 


945 



Now we provide a general method for carrying out summation over a's. We 
express the sum from 1 to p over each of the n a's as a sum over the partitions of 
n. In our example of the four-point function, we have reduced the calculation of, 
for example, the numerical coefficient of the middle term in Eq. (6), to the sum 

p 

Yl (2<^™-l) (22) 

Instead of summing each a from 1 to p, we can break up the sums into a sum over 
single values of each a, for example 

E = E + E+--- + E (23) 

ai=l ai=l ai=2 cei=P 

When we carry out multiple sums, we have to take account whether there is a {26—1) 
factor connecting a pair of a's and which sets of a's match and which don't match. 
If there is no {26 — 1) factor connecting a pair of a's there is a factor 1; if there is 
a {26 — 1) factor connecting a pair of a's there is a factor —1. Next, we select any 
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specific term, and sum over the set of, say, h a's, that do not appear in the summand. 
Obviously summation over these /i a's leads to a multiplicative factor of in the 
coefficients. Now we are left with the summation over the remaining {n — h) = f 
a's. To make further progress, we begin with the case oi h — corresponding to 
the maximum number of crossings. Here not only all a's appear in the summand, 
but each one appears for the equal number of times. This symmetric distribution 
of a's in each summand leads to considerable simplification. Our further discussion 
applies only to this "maximally saturated" case, although the use of partitions with 
restrictions may apply generally. 

For this symmetric case wc group the numbers from 1 to n into groups that 
arc equal and unequal. This leads to unrestricted partitions of the number n. Let 
the p{n) partitions of n be labeled X^^\ X^'^\ • • • , A*^^*^"''^. For a general A*^*^ let 

A(^) = (AS^^..Ai^)A^^)...Ar...A;^)...A;^)...A!^)...Af)) (24) 

where m[^^ is the number of times the part A^*^ of the partition A^*^ occurs, k^^^ — 
YliJi "^1*^ is the number of parts of the partition A^*^ l'^^^ is the number of distinct 
parts, and the weight is |A^*^| = Z^fjii) — ^- Then 

-1 n[(Af)!)™'^ mf)!] 

i=l 

(25) 

where is the total number of 'n-tuples' {ai,a2, ■■■,(y.n) when each a varies from 
1 to p. In our case, the /(ai, • • • ,an)\\(s) is the product of the factors {2S — 1) 
associated with the crossings discussed above. This product is either 1 or -1 for a 
given partition. For the case of no crossings, whose n-vertex graphs consist of n 
disconnected points, f{ai, ■ • • , an)\x(s) = 1 and we find 

i=l 
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as we expect. The summand for any 'n-tuple' of (ai, 02, .., On) always reduces to 
either 1 or —1. The partitions introduced above determine this sign. As noted 
above {2Sai,aj — 1) = 1 if = and —1 otherwise. If in any given partition A*^*-*, 
ai is not equal to aj for A^'^^ number of times, the summand would be 
This is true for all the possible groupings of a's according to the partition A^^'*''. 
This is a consequence of the fact that in the most saturated case the a's have a 



symmetric distribution in the summand. From Eq. (25), we get the total number 
of configurations E^'^^ associated with the A^*-* as 

^ . x(') = — (27) 



i=l 



The summand contains a product of n{n — l)/2 (25aj, aj — 1) terms, and in 
the partition A'-^^ the number of times for which equals is 

Ms)^Y^^) \ -1) . ;^(^)>2. (28) 

Consequently, 



i=l 



and the coefficient 



p(n) 

E n (25..., - 1) = (30) 

ii,---,"n i<j;2<{i+l j}<n s=l 



As an illustration, we evaluate the coefficient C|| for the most saturated graph 
corresponding to = 12, viz., 

= E n (2^".-.-i) (31) 

01,.. .,.6 i<j;2<{i+l,j}<6 

It contains a product of 15 {26 — 1) factors or a sum of 2^^ terms. The relevant 
parameters for calculating C|| are given in the TABLE II. Using Eq. (32) and 
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Table II, we get 

CI5 = -p{p - 2)2 - 26p2 + 152p - 
For bosons, p = 1, and C|| should, and does, reduce to 1. 
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TABLE II 



s 
















1 


[6] 


6/1 


- 


- 


1 


6!p!/{6!(p-l)!} 





2 


[5,1] 


5/1 


1/1 


- 


2 


6!p!/{5!(p-2)!} 


5 


3 


[4,2] 


4/1 


2/1 


- 


2 


6!p!/{4!2!(p-2)!} 


8 


4 


[4,1,1] 


4/1 


1/2 


- 


3 


6!p!/{4!2!(p-3)!} 


9 


5 


[3,3] 


3/2 


- 


- 


2 


6!p!/{3!22!(p-2)!} 


9 


6 


[3,2,1] 


3/1 


2/1 


1/1 


3 


6!p!/{3!2!(p-3)!} 


11 


7 


[3,1,1,1] 


3/1 


1/3 




4 


6!p!/{3!2(p-4)!} 


12 


8 


[2,2,2] 


2/3 






3 


6!p!/{2!33!(p-3)!} 


12 


9 


[2,2,1,1] 


2/2 


1/2 




4 


6!p!/{2!^(p-4)!} 


13 


10 


[2,1,1,1,1] 


2/1 


1/4 




5 


6!p!/{2!4!(p-5)!} 


14 


11 


[1,1,1,1,1,1] 


1/6 






6 


6!p!/{l!66!(p-6)!} 
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3 Summary and Conclusions 

Wc used Green's ansatz to evaluate the vacuum matrix element of an arbitrary 
product of parabose fields of order p. We gave a graphical method to reduce an 
N — 2n point function to a sum over products of n 2-point functions. We expressed 
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the coefficients of eacli tliese n 2-point functions as a sum over n a^'s (1 < ccj < p). 
The corresponding graph uniquely determines both the n 2-point functions and the 
coefficients that occur in the sum. 

Subsequently, we gave a combinatorial method, based on the partitions of the 
number n, to determine the coefficients. We illustrated this method by the exphcit 

determination of the coefficient for most saturated graph corresponding to the 
N = 12 case which involves a sum of 2^^ terms. 

4 Figure Captions 

Figure 1. Case (1) of a 5-crossings graph for an 8-point function. The algebraic 
expression associated with this graph is given in Eq. (20) 

Figure 2. Case (2) of a 5-crossings graph for an 8-point function. The algebraic 
expression associated with this graph is given in Eq. (21) 

Figure 3. Case (3) of a of 5-crossings graph for an 8-point function. The associated 
algebraic expression for this graph is given in Eq. (22) 

Figure 4. Case(4) of a of 5-crossings graph for an 8-point function. The associated 
algebraic expression for this graph is given in Eq. (23) 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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